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F* (§ ) 
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J 
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see p. 39 
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n ’ ’ 
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y + 1 M 2 * 

" 2 b 2 


k 

see p. 41 


L 

- 2 (K + 2 

M 2 - f ) 



V 


V n * 

0 . 1 

see ( 5 . 18), (D. 8 ) and (D.16) 

M 


free -stream Mach number 

m l 


local Mach number 

— * 
n 


unit vector normal to the shock 

orxi 


see p. 10 

3? 

3 

II 

1,2,3. 

arbitrary real numbers 

P 


row vector (p 1 , p 2 , p^) 

Q 


see p. kO 

q 


flow velocity vector 

q i 


= u 

q 2 


= V 

q 3 


= rV 

q 


Bernoulli constant (see p. 5 ) 

R(x) 

, 

r = R(x) ss eR(x) is the body contour 

> 

R(x) 

J 

axisymmetric case 

R (x ) 
o' o' 


r Q = R q (x o ) is the image of the body 
in the x r -plane 

r 


cylindrical coordinate 

a 

a 

0 , 1 , 2 , . . . 

see (3.3), (5.29) and (5.34) - (5.35. 

r oB 


see subscript B below 

r 

m 


see Fig. 6.2 

U 


q = Ue + Ve + We 0 
x r 9 

U 

00 


free-stream speed 

if 

3 

II 

1 , 2 ,... 

see p. 19 and (5«l) 

u 


U = U (1 + u) 
00 . 

u 


row vector (u,v,w) 

u o 


see ( 3 . 21 ) and ( 5 - 29 ) 

u n , n = 

1,2,3,... 

see ( 3 . 3 ) 

“b 


= u qB see subscript B below 



V n , n — 1,2,3,... 


v n y n — 1- j 2 ,3 » • 


B 


W 


W n t n _ 1 , 2 , 3 , • • * 


w 


x n » n — 0, 1 ,2,3 , . , 


oB 


x (r) 
s v ’ 


n 


n = 1,2,3. 


yi, 

see U , • 

see p. 19 and ( 5 . 1 ).. 

V = U v . 

CO ' * 

see (3 .22.) and ( 5 . 29 ) 
see (3.3) , , • 


•=V Q g see subscript B below 


see U 

see ( 5 . 1 ) 5 

* v~W = U w 1 ■ ■■■" 

oo - 

see ( 5 . 29 ) 

cylindrical coordinate 
row vector (x,r,9) 

see (3. 3), ( 5 . 29 ) and (5.34) - (5.35); 
in Appendix B x ^ = x, = r and x^ 

see subscript B below • >••• ' - - 

• x- = x (f ) is a,sh'ock surface 

s' . 

see. p. 41 

row vector (y r , ^ 2 , v y^) ' * 


Greek alphabet 


O' 

V 

e 

C 


T) 

e 


® n » n — 0, 1,2,3,... 


sin a = 


•i j 

_ 1 _ 


= c /c 
p' y 


overall thickness ratio of the body 
characteristic variable (see p. 44 with 

y-i = 5 » y 2 ' = H 311(3 y 3 = C ) 

characteristic variable (see £ and p. 
cylindrical coordinate 
see ( 5 . 1 ) 
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flow deflection angle (V = U tg ,?) 

*1 


see ( 5 * 19 ) 

H ? 


see (C.7) 

x(s) 


see (3.51 ) 

V n = 

1 , 2 ,,,#, 6, 

see (3. 8), (3.9), (3. 15), (3. 16), (3. 36), (3. 37) 
(3.40) and (3.53) - (3.55) 

X (n) ,n = 

1,2,3. 

see p. 45 ' 

2 v- 1 

M< 


ll 

11 

5 


characteristic variable (see £ and p. 38) 

a 

i 


T 

7 

see Fig. 6.1 

4> 


velocity potential (*} = V<f) 

—* -4 — • 

0 


ue + ve + we A » = Vj 
x r 0; - , 

0 


see ( 3 . 25 ) and ( 5 . 29 ) 

V n = 

1,2,3,... 

see (3. 12) . and (3.3 1) 



= 0 n see subscript B below , ; 

X 


= x - Br ; . . , 

x n > n = 

0,1,2,... 

see p. 19 and ( 5 . 1 ) 

,( r ) 


x = \|i(r) is a shock surface 

,(*) 


\|r(x) = const, is a characteristic surface 
see (b.9) 



see (B. 1 1 ) 
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o 
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Subscript B 

3_ 

3x 

a_ 

dy 

a_ 

9p 

a(x) 

a (y) ’3 (x 1 

a (x , r ,e ) 
a(s.Ti.c) 


0) 

( 2 ) 

(a) 


» y 2» y 3 j 

» x 2 jX^j ) 


denotes quantity in front of the shock 
behind 

the corresponding axisymmetrio expression 
values taken at the body contour (see p. 15 ') 

-3 3 3 


1 a a 0 \ 

row vector 


/3 a a \ 

Wi ’ Sy 2 ’ 3y 3 / 


I o o 0 \ 

\3p 1 ’ 3p 2 ’ 3pJ 


and 


> are functional determinants 
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1 . INTRODUCTION 


In the well known Whitham theory [8] for predicting the strength of 
the sonic boom at large distances from a supersonic vehicle only 
such nonlinear effects are included which produce a distortion of 
the signal as it travels away from the vehicle. These have a cumula- 
tive effect at large distances and are therefore important even for 
very small disturbance levels. Whitham's simple first-order rule, 
namely that the perturbation velocity field is the same as the one 
obtained by linear theory but displaced in the streamwise direction 
so as to take into account the distortion of the downstream Mach 
cones, will usually be satisfactory for the very slender configura- 
tions of interest for commercial SST's when flying at moderate 
supersonic Mach numbers. For not so slender configurations or high 
Mach numbers (above 3, say) higher-order nonlinear effects can be 
expected to become important in regions close to the vehicle. For 
larger distances, say a few body lengths away, linear acoustics 
would still be expected to hold, although the signal strength would 
not be correctly given by the first-order theory relating signal 
strength to vehicle geometry. 


To assess the importance of such higher-order nonlinear effects in 
the near field, a second-order theory for the complete flow field 
would be of great value as it would allow one to estimate the errors 
inherent in the first-order theory. Such a theory for the near field 
of a body of revolution at zero angle of attack was obtained some 
time ago by Van Dyke [7] by the use of a regular perturbation method. 
This solution is not uniformly valid in the far field domain — or in 
regions influenced by discontinuities in the body surface slope — 
since the second-order solution thus obtained has the same set of 
characteristics (free-stream characteristics), and hence influence 
regions, as the first-order (linearized) solution. Therefore, 

Van Dyke's theory needs to be modified so as to account for the 
change in the characteristics. A systematic expansion procedure 
yielding a uniformly valid solution in the whole flow field has been 
proposed by Lin [2] and further extended by Oswatitsch [3]. In this, 
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both families of characteristics are readjusted at' each step. In 
an earlier publication* [ 13 ] , it was shown that this procedure could 
be used to cast Van Dyke's second-order solution for axisyrrenetric 
flow in a‘ uniformly valid form for moderate and large radii. The 
present report shows how the uniformly valid solution can be obtained 
directly by aid of the Lin-Oswatitsch formal expansion procedure. 

Van' Dyke 's results are obtained when the solution is expanded in a 
Taylor series around the free stream characteristics. Thus, an 
extension of Whitham's Q8] rule to second order is presented. 

The hew theory is also' modified to allow for discontinuities in body 
surface- slope.' For the determination of shock wave location a simple 
extension of Whitham's "rule of equal areas" is presented. 

The three-dimensional far field is treated in a similar manner using 

1 a perturbed cylindrical Coordinate system. For' the three-dimensional 
mid and near fields, the use of characteristic variables does not 
appear to lead to a simple second-order solution. The far-field 
solution is used to generate formulas relating flow deflections to 
the Whitham F-function. With the aid of such formulas a new experi- 
mental technique was devised [14] for the determination of the sonic 
boom strength from wind tunnel measurements which, particularly for 
high Mach numbers, can be carried out fairly close to the body where 
the disturbance levels are large enough to give good experimental 

-accuracy. 

.2. GENERAL ASSUMPTIONS, BASIC DIFFERENTIAL EQUATIONS AND 
CHARACTERISTIC COORDINATES 

We study the steady supersonic flow field around a fixed rigid body, 
placed in a homogeneous parallel stream of a polytropic ideal gas. 
Viscous forces, gravity and heat conduction are neglected. Cylind- 
rical coordinates x,r and 9 are introduced such tha-t the free 
stream is in the positive x-direction and the body is situated in 
the half -space x 2: 0. Further, it is assumed that all streamlines 
start from x = -CD. 
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The entropy .change across a weak shock waye. -is of the order -of the 

shock strength to. the third power, i.e. of the order of the cube 

of the velocity perturbation jump across the wave. In the second- 

order theory considered here only, contributions from. terms. which. _ 

are quadratic in the velocity perturbations, are retained , . except. .. . 

near a slender body of revolution, where some, terms which are 

formally cubic in the v-perturbation. ; are included ; For a smopth ; • 

pointed body, however, the bow shock strength vill.be of the order. 

2 

u also near the body, so that entropy changes .may .still be neg- 
lected in this region. If the body has slope discontinuities of 
order e, where e is the overall thickness ratio of -the body,., r- 
the perturbation velocities and hence the shock strength will be 

of order e, so that entropy changes and . the entropy, gradients. 

' : 3 

become of order e , which may be neglected in the second-order 

2 

theory then retaining only terms of order e , . Thus.,-, to the approxi- 
mation considered, we may neglect the. entropy, changes and .hence con- 
sider the flow to be isentropic and irrotational throughout the flow 
field.. . We may thus startfrom the equations, for irrotational flow . 



with the velocity of sound c given by 

^ . * .'«*) . - 

(Bernoulli's equation). 
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Here q is a constant, y = C^/C^ and U, V and W are defined by 


q = U e x + V e r + ¥ e Q 


q being the flow -velocity vector. 


The above equations together with the shock relations and the con- 
ditions for the flow to be tangent to the body surface and to tend 
to the free-stream flow for x - co specify the flow iri the whole 
field. 


In the following we consider x and r as dimensionless coordinates, 
obtained through division by a suitable reference length such as the 
overall streamwise dimension of the body. This does not change the 
form of (2.1) - (2.4). We shall also use dimensionless disturbance 
velocities u, v and w defined by 

U = U (1 + u) 1 

V = U v 

CO 

V = U w 

00 J 

U denoting the free-stream speed. 

00 

In the axisymmetric case (with the x-axis as symmetry axis) we have 
9/^9 = W = 0 so that (2.3) arid (2.4) need not be considered. In 
that case it is advantageous ( [2] , £3] ) to introduce characteristic 
coordinates | and T| as new independent variables (see Appendix A). 
If we choose § and T] such that | = const, is a downstream cha- 
racteristic and T] = const, an upstream characteristic the trans- 
formation ( 5 » T| ) -* (x,r) can be written 

O+u-Bjv)!! + (b l+ E l u + v)|| =0 (2.5) 

k ( 1+U+B L V ^ - K +b l u - = 


0 


(2.6) 
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and equations (2.1) and (2.2) are then transformed into 
( V Bjj-i - v )|^ + ( 1 + u + B^ )|-^ + ( 1 + u + B l v ) ( 1 + GjV )— ~ ■ = 0 (2.7) 


(Bl+ Bl u +v)|^ - ( 1 + u - BjV )|^ (1 + 

Here Bj^ = ^M^- 1 and G L = M 2 (l - :1 ~ 
Ihdex L denotes local values and M 


u -B l v)(1 +G l v 2 )J|^= 0 (2.8) 

M 2 (2u + u 2 )- M 2 V 2 ) . 

the free stream Mach number. 


From Bernoulli's equation we obtain 


B l = B 


-1 


with K = \ 

B 


1 + (K- M 2 )(2u + u 2 + v 2 )( 1 - 4^1 M 2 [2u +u 2 + v 2 ] ) ] 

.4 


and B = VM 2 -1 . 

Each of the equations (2.5) - (2.8) contains derivatives with respect 
to only one variable. This simplifies the integration of the dif- 
ferential equations considerably. 


In the general three-dimensional case one speaks of characteristic 
surfaces and bicharacteristics belonging to the system (2.1) - (2.4) 

(see Appendix B) . For the special case of undisturbed flow, i.e. 
q = e x — which is a trivial solution of (2.1) - (2.4) — the down- 
stream Mach cone x - Br = const, and the upstream Mach cone x + Br = 
const, are characteristic surfaces and the meridian plane 9 = const, 
intersects the downstream Mach cone along a bicharactoristic (see Appendix B) 
In our study of the general three-dimensional case — with the only restric- 
tion that the disturbance caused by the body is small — it appears to be 
advantageous to introduce new independent variables §,T| and Q such that 

5 = const, is a characteristic surface and a disturbed down- 
stream Mach cone (x - Br« const.), 

T) = const, is a characteristic surface and a disturbed up- 
stream Mach cone (x + Br » const.), 
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Q = const, intersects 5 - const along a bicharacteristic 
and is a disturbed meridian plane ,(9 ** const.)'. 

In the following, the variables § , T) and Q' will be referred to 
as characteristic coordinates. The transformation (^jTI.C) - * (x»r,9) 
can be written formally (see Appendix B) 


( 



9 x 

3.S 



M, y ) 


0 


(2.9) 


<■ f 2 ( r ’ u > 371 > 5 M* v ) “ 0 

v ^ 3 ( r ’ u ’ §tT’ Tc 5 M> v ^ = 0 


( 2 . 10 ) 


(2.11) 


Also the equations (2. l) - (2.4) can be transformed to equations of 
the form - ; ’ - 


I f^(-r,u. 


9x 

a?’ 


dx 9x 9u 9u <. 

aV at* a?’ 3C 5 M ’ Y } 



( 2 . 12 ) 


/ - 9x_ 9x 9x 9u. 9u 

5 lr,U ’ 9?’ 311' 9C* 311’ 3£ 


M, Y ) 


0 


(2.13) 


/ 3x 9x 9x 9u 9u 9u \ . n 

o' r » 35» 91] ’ 9C’ 3?’ 9 T| * . 3 C ' 


(2.14) 


V 


. 9x ,9x 9x du.au. 9u \ >, 

7 (r,W ’ 9|’ 911’ a C * 9 5 ’ 9 "H ’ 9C ; . 


(2.15) 
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Here x and . u denote the vector quantities (x,r,9) and (u,v,w), 
respectively,, and the f^'s are complicated expressions given by the 
left-hand sides of equations (B.23) - (B.29). The following can be 
said about (2.9) -(2.15)5 

Equations (2.9) and (2. 10) are the conditions for. T) = const, and 
5 = const., respectively, to be characteristic surfaces. In case of 
axisymmetry they are equivalent to ( 2 . 5 ) and (2.6), respectively. 

Equation (2. 11) is the condition for £ = const, to intersect £ = const, 
along a bicharacteristic. In the case of axisymmetry equations (2.12) 
and (2.13) are equivalent to (2.7) and (2.8), respectively. Similarly 

(2. 14) and ( 2 .I 5 ) are just (2.3) and (2.4), respectively, (f^ does 
not contain 3v/^£, 3 v/3t], and dv/3£ and f^, does not contain 

d u/^ 5 , d u/& 1 ) , and d u/£) £ ) . 

As in the axisymmetric case, absence of certain derivatives in (2.9) - 

( 2 . 15 ) helps simplify the integration procedure. 

3. A UNIFORMLY VALID SECOND-ORDER SOLUTION FOR AXISYMMETRIC FLOW 
AROUND A SLENDER BODY OF REVOLUTION 

Let the body considered be located in the half-space x S 0 with- 
the positive x-axis as axis of symmetry and let its contour be de- 
scribed by 

r = R(x) = eR(x) , : ( 3 . 1 ) 

assumed to be pointed, and is allowed 
of order e.’ The flow must satisfy 
body surface, , 

v = R ( 1 + u) (3.2) 

and must be undisturbed upstream of the body. In addition,, shock 
waves give rise to boundary conditions on surfaces x = x g (r) which 


where R = 0(l). The body is 
to have slope discontinuities 
the tangency condition at the 



cannot be specified a priori. It turns out that shocks lead to a 
nonuniqueness in the transformation (§ , T| ) -» (x , r ) which requires 
special treatment as shown in the subsequent section. 

Jn attacking the problem for axisymmetric flow we shall first de- 
termine a second-order solution valid in the mid-field region 
#j 

r ord (l) and then extend this to regions near the body and far 
from the body. 


For the mid-field expansion we assume, following the procedure 
proposed by Lin [2] and further extended by Oswatitsch [3] , 


u 

= u i 

+ 

u 2 

*4* « • • 

V 

= v i 

+ 

v 2 

+ • • • 

X 

n 

X 

0 

+ 

X 1 

+ • • • 

r 

0 

u 

n 

+ 

r i 

+ • • • 


(3.3) 


where the different terms are to be regarded as functions of the 

characteristic coordinates 5 » T| • Here, x_ and r_ are of order 

2 


and 


unity, u i , v 1 , x., 

(3.3) into (2.5)-(2.8) we obtain 


1 


of order e 


o o 

Upon substitution of 


rr (x + B r ) = 

v o o' 


o 


(3.4) 


<L_ 

an 




o 


(3.3) 


a? 


(B u 1 




a? 


o 


(3.6) 


a_ 

an 


(B u 1 



ii ifo 

r o ^ 


0 



r 


ord (l) 


we mean 


{ 


r 

-1 

r 


0 ( 1 ) 


(3.7) 



The solutions of (3.^) and (3*5) read 
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x 

o 


+ B r 

o 


X (Tl) 


(3.8) 


X 

o 


B r 

o 




(3.9) 


Without loss of generality we may choose X. -j = T| and \g = % • The 
characteristics will appear as straight lines in the x. :r o -p.lane , 
and it may be convenient, following Oswatitsch [3], to use 


X o = i (Tl + S) 


(3.10) 


r o = ' 2 V (71 ~ «> (3 * 11} 

in lieu of 5 and T) as characteristic coordinates. 

Substituting (3.10) and (3.1l) into (3-6) and (3.7) we obtain the 
usual linearized equation for u^ and , i.e. 


2 

0 d 0 -1 1 S0 1 

- B + - 

dx 00 


aV 


dr 


= 0 


(3.12) 


where u 1 = d 0 1 /^x Q , v. ] = d 0 1 /^r Q . 


We proceed now to next order. From 

a sXq 

af ( x i +B r i) + ( u i“ B v i)al“ + 


^(xi'Br,) •+ ( u i +B v i)arT " 


( 2 . 5 ) and ( 2 . 6 ) we obtain 
P d r 

(B[K - B ]u 1 + V.,)— = 0 


(B[K-B 2 ]u 1 - Vl )^ = 


(3.13) 


0 (3.14) 
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It, is readily verified using (3.10) 
(3.13) and (3.14) can be written 

x i + B r.j + M 0 1 + K r Q ( B u 1 + v ^ ) 
B r 1 + M 2 0 1 - K r Q (B u 1 - 


- ( 3 . 12 ) that the solutions of 


x 3 01) 

(3.15) 

X 4 (?) 

( 3 . 16 ) 


The discovery of the explicit solution first given in [ 12] makes pos- 
sible considerable simplifications of the results. 


Becuase of the condition that disturbances vanish upstream of the 

body so that the Mach lines then become those of the undisturbed 

f . •' - • . • ' - 

flow we select X 3 5 0. Also, the solutions for u ^ and v^ turn 
out to be particularly simple with the choice \ ^ h 0. Thus, we have 


x, = -Kr v,(x,r)-M0 1 (x,r) 
1 o 1 ' o o' 1 O o' 


(3.17) 


’1 = - K r o U 1 ( x o’ r o> 


(3.18) 


Now we return to ( 2 . 7 ) and (2.8), Substitution of (3.17) and (3.18) 
and retention of the next order terms yields 


[If < B u 2 + v 2 ) + ~r 1 + L b(k - b2)u i _ v il aT 


^ V 1 V 1 9r o ^ U 1 

+ ( U 1 + B Vl )^~ + (u, + B Vl )— gy - K Vl 


5 

dv 

H 


= 0 


( 3 . 19 ) 


h. ( B 


.) - 


- ( u i “ B 


12 Ho 

r o ^ 


] + [b(K-B 2 )u 1+Vi ] ^ - 


9v 1 V 1 3r o 

v-^“ - ( U 1 -B Vl )— r=- + K v 


^1 
1 aii 


= 0 


(3.20) 


With the aid of (3.6) and (3.7) and with introduction of the new 
dependent variables 
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(3.21) 

( 3 . 22 ) 

we find that ( 3 . 19 ) and (3.20) may be written in the following simple 
way : 


• ' „ 2 - 2 

u = u, + u„ - M u, 
o id \ 

v o = V 1 + v 2 -<^ 2 + K) u 1 


d_ 

a? 


(B u Q 


a_ 

*t\ 


(B u 
x o 




+ 



v dr 
o _o 

r o . 


= o 



(3.23) 


(3.24) 


Thus, a comparison with (3.6) and (3.7) shows that u q and v q are 
obtained from the solution of the . linearized equation 





0 


(3.25) 


with u = b<b /bx , v = b<t> /b r . Therefore, the second-order 
o o' o o *o' o ’ 

solution can be expressed directly in terms of the solution of the 


linear equation ( 3 . 25 ) by 

u = (1 + M 2 iL) u Q (3.26) 

v = ( 1 + [K + M 2 ]u q ) v q (3.27) 

where u and v are considered functions of x and r , and 
00 00’ 

x and r are given, by 

x = x q - M 2 0 o - K r Q v q (3.28) 

r = r Q (1 - K u Q ) •- • 1 • (3.29) 


One can also readily show that the velocity potential to second 
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order is given by 


0 = 0 -K r u v 

o o o o 


(3.30) 


Next, we shall extend this solution to the region close to the body, 
r = 0(e), by making use of Van Dyke's [7] solution. Through a regular 
expansion procedure Van Dyke found the solution for the velocity 
potential 0 = 0 1 + 0^ to be 


2 M 2 

0 2 = u.j (m 0 1 + Krv^) - Tj-rv., 


(3.31) 


from which the following velocity components can be obtained: 


u = u 1 + u 2 


2 9u 1 du „ 

= Ul (l + M u 1 + Kr r— )+ 5 ^(M 0,, + Krv,) - 


3 m 2 9U 1 2 

i M r V 1 


(3.32) 


2 9v 1 dv 2 

v = v i + v 2 = Vl [ 1 + (M + K) u 1+ Kr§^-] + ^-(M 0 1 + Kr Vl ) - 


M 2 / 3 o 2 &V b 

” TT ( v i ° rv i 5^-) 


(3.33) 


These are to be compared to the velocities, obtained by expanding the 
uniformly valid solution for r = ord ( 1 ) obtained above. 

Su 1 au 1 

Setting u,(x , r ) = u, + (x - x )- + (r - r )^ + , r = r + . . 

1' o’ o' 1 ' o/3x ' o'3r o 

etc,, we obtain 


2 3u 1 9u 1 2 

= Ul (l + M u, +Kr + + Krv i) 


(3.34) 


2 9v 1 dv. 

v = Vl [ 1 + (M + K) Ul + Kr ^-] + j^-(M 0 1 + Krv,) 


(3.35) 


which differ from ( 3 . 32 ), (3.33) only in that the triple product terms 
are missing. Van Dyke £7] pointed out that these will be of the same 
order as the quadratic ones near a slender body for r = 0(e). For 
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= ord (l), however, they are negligible and of the same order as 
contributions from other terms already neglected in the differential 
equation. It is possible to make the solution (3.26) - (3.29) uni- 
formly valid also for r = 0(e) by modifying it as follows: 


* - *„ - K r o v o - <^ Br o v o 2 + l x 5 <?) (3 - 36) 

r , r c (1 - K u c ) - l M 2 r o v/ - jg XjfS) (3.37) 


V 


[1 + (K + M) u o 



2 , 

V V 

o J 


o 


(3.38) 


The expression for u q is the same as before. The additional terms 

involving the as yet undetermined function \ _ , and the one with 

2 ^ 
v q in (3.36), are everywhere of higher order and incorporated to 

avoid discontinuous behaviour of x when u and v have dis- 

o o 

continuities at a body slope discontinuity. The term r has been 
included in order to eliminate the spatial derivative terms in the 
triple product contribution to u and v given by the last terms 
in (3.32) and (3.33). It becomes comparable to the previously 
derived term -K r u only for r = 0(e) and is therefore at most 

3 ° ° 

of order e . In constructing the new term we have freely taken 
advantage of the s lender-body behaviour of v q . Thus, for example, 
for r = 0(e ) , dv./fer ~ dv /dr c* -v / r 

To determine > let the body contour in the transformed plane be 

defined by 

r = R (x ) 
o o o' 

By choosing 

X 5 = -2BK^( S ) R 0 (x oB ) - | M 2 B v l R o (x oB ) Ci.hO) 


(3,39) 


oB’ r oB 


where 


"B = U o (x oB’ r oB^ 


V B = V o( X oB’ r oB> X 


are the 



14 

coordinates for the intersection of the particular characteristic 
5 = const, with the body contour (see figure 3.1) we can ascertain 
that both x q and r Q become continuous at points of body slope 
discontinuities. ' 

r / 



Figure 3.1, 

This follows, because in the immediate neighbourhood behind the 
body comer the flow must behave locally like a two-dimensional one 
with the velocity component jumps to lowest order satisfying AUg = 
- AVg/B = AR ' /B , where AR' is the slope discontinuity. From 

the' expressions (3.36) and (3.37) it then follows that a continuous 
body contour in the x q , r Q plane leads to a continuous' contour in 
the physical plane with 


r = R = R [x ^(xj] 
o L oB v 


(3.41) 


and x given implicitly by 


x ..= X oB - M 0 B - KR o( v B + Bu B ) 


(3.42) 


(Continuity in x ' follows since the combination v q + Bu q is con- 
tinuous, as well as 0 g .) Thus, a solution, uniformly valid to second 
order everywhere for r — 0( 1 ) , is given by 


u = u ( 1 + M u ) 
o v o' 

= v[1 + (K+M 2 ) u o - v 0 2 ] 


(3.43) 

(3.44) 


V 
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2 

x = x o - K r o v o-^o* BKR o “b'T - ( r o v o 2 - R o V B 2 ) (3.45) 

r = r o (l-Ku o ) +KR q ^-^(r, v q 2 - R q v 0 2 ) (3.46) 


A comparison with the previous result shows that the added terms 

give contributions which are of hi^ier order everywhere with the 

exception of the added term in v , which gives a contribution of 

the same order near the body as the term (K-M )u q . For r = ord (l) 

it becomes of higher order and may be neglected. The additional 

3 

terms in the expressions for x and r are of order e everywhere 
and are only important for representing the details of a non-smooth 
body in the transformed plane. It should be noted, that for a body 
contour with slope discontinuities, the transformations x = x(x Q ,r o ) 
and r(x Q ,r o ) as given by (3.45) and (3.46) will be discontinuous 
outside the body along the downstream characteristic from the corner. 
For a convex comer (expansion) there is a negative jump in v q and 
a corresponding positive one in u Q giving a positive jump in x 
and a negative one in r as illustrated schematically in figure 3.2 a). 





a) Convex comer ■ b) Concave comer 

( expansion) ( compression) 


Figure 3.2. Coordinate transformation near discontinuities. 
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A double point A q A^, on the upstream and downstream sides, respect- 
ively, of the characteristic 5 = const, from the corner thus becomes 
transformed into a line along T] = const, as shown. The gap between 
A and A' in the physical plane must of course be filled with an 
expansion fan, as is well known from the two-dimensional theory. 

In the present problem the actual variation of flow quantities across 
the fan can be determined by considering the sharp comer as the 
limiting case of a rounded one as the radius of curvature tends to 
zero. For a concave comer, on the other hand, the downstream point 
A ' is transformed to an upstream one in the physical plane (see 
figure 3-2 b)) and a problem of non-uniqueness thus arises, since one 
and the same region in the physical plane might correspond to two 
regions in the transformed one. This difficulty is to be resolved 
by the insertion of a shock wave to remove the overlap region in 
the transformed plane. Such regions may of course also appear even 
for a body without sharp concave comers, and the problem of deter- 
mining shock wave locations will therefore require a more general 
treatment given in the following section. 

With the uniformly valid solution near the body thus derived we may 
proceed to consider the boundary condition on the body surface. By 
introducing the second-order solution (3-43) - (3.46) into the tangency 
condition (3-2) we find upon neglecting higher-order terms, 

v q [ 1 + (K + M 2 )u q - \ v q 2 ] = (l+u o )R'(x) (3.47) 


which, together with the initial condition of zero disturbances up- 
stream of the body, provides the boundary condition for the linear 
equation ( 3 . 25 ). As the transformed body contour is continuous, 
the problem of finding the transformed body contour in the x Q ,r - 
plane can be solved by iteration starting from the ordinary first- 
order solution with the linearized boundary condition 

Vl (x ,R(x)) = R'(x) 


(3.48) 
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It remains to be studied the behaviour of the solution in the far 
field, r = ord (e^*), § = 0 (l). In Appendix D it is shown that by 
assuming 

4 tt 8 

u = e U^ + e Ug + . . . 

4 __ 8 

v = e + e Vg + • • . 

x - Br = X = X Q + ^ X-, + ••• 

4 4 

e * 2Br = H = H q + e H^+... 

where U^, V^, and are functions of § and T] , all of 

order unity, and H q = ord (l), one finds that the solution may be 
written 


+ T 7^375 ^<5')d5W(M 2 -|)^ + 0( e ' 2 ) 0.49) 


(2Br) 


v = - Bu + 


B 


(2Br) 


J/2 


2 

^ P(§')d5' " BK fi^ + °( e12 ) (3-50) 


2 K 


x - Br = § - K \J 2Br' F(| ) + 


M - 


/2 B? 


^ F(5')<*5' 


(3.51) 


+ LF 2 (|) In Br + \(s) + 0(e 8 ) 


^ T K lv 3 „2 2K x 

where L = -g- (K + M - — ) 

M 


From the mid-field solution expanded in the same manner one finds, 
using the far-field expansions for u q , v q and 0 Q given in 


the 
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o 

Appendix C, and noting that from (3.37) r 0 — r + Kru Q + 0(e J ) , 
same expressions for u and v as above, but that ( 3 - 51 ) is re- 
placed by _ 

m 2 k * 

I — * M - r f 

x - Br = 5 - K V 2Br F(? ) + — \ F(^ « )d§ ' + 

/2Br ' ^ 

o 

+ |m 2 b 2 f 2 (|) n 5 ( 3 . 52 ) 

The last two terms may be absorbed in the function X(?) of ( 3 . 5 1 )- 

However, the term proportional to In Br which is the result of 

higher-order cumulative effects in the characteristics in (3.5 1 ) 

has no correspondence in (3.52). We thus need to complement the 

2 

midfield solution with terms which behave like LF In Br in the 
far field, but which will only give higher-order contributions near 
the body. The added terms, however, must not lead to discontinuities 
in the transform of the body contour. Hiis can be achieved by re- 
placing (3.36) and (3-37) by 


x 


X 


o 


-Kr v -M 2 0 + jr M 2 B r v 2 + 

o o o 4 o o 

+ LBr (in B r )u 2 + ■=■ X , 
o'- o' o 2 6 


(3.53) 


r (1 
o' 


K u 


- i 


M 


- L 


r ln(Br )u 
o ' o' o 


2B X 6 


(3.54) 


where 


X* = - 2BK u n (|) R (x J - | B v 2 R (x J - 
, 6 o' oB' 2 B o' oB 7 

- 2 L B R q m(B Rji^ 2 (3.55) 

The transformation (3.53) - (3.55) together with (3.43) and (3.44) 
provides a solution which is uniformly valid to second order every- 
where. The new terms in (3.53) and (3.54) give contributions to x 

4 

and r which are at most of order e . In practice they are likely 
to be unimportant since they lead to contributions in the mid field 

g 

of order e in the far field. Finally, it should also be remarked 
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that a uniformly valid solution can be constructed in many different 
ways and the one given here was selected with computational simpli- 
c ity in mind . 


4. DETERMINATION OF SHOCK WAVES FOR AXISYMMETRIC FLOW 

At a shock wave, the following conditions hold (see, e.g., [ 1] , 

P. 274): 

5<')x3< 2 ) = ♦ (l-u 2 )(5 (,) *nf 

q(')x.s = 5 |2 ) x: 

where and q^ 2 ^ are the velocity vectors ahead and behind the 

shock, respectively, p, 2 = (y-l)/(Y+ 1 )» q 2 = U 2 + V 2 + W 2 + 2c 2 /( y -1), 
c the local velocity of sound, and n is the unit normal to the shock. 
Let the shock in the axisymmetric case be defined by x = <|i(r). Then 
(4.1) gives 

(S) 2 = B 2 ^[ 1 *(i) 2 j(u('),uW,uO) U ( 2 ), 

+ v^M 2 )- 2M 2 g + v^-M 2 ^ 1 * g+v* 1 *) (4.3) 

Expanding this result for small values of the perturbat ion velocity 
components and retaining only terms of zeroth and first order in. 
consistency with the second-order theory for the mid-field, we find 

g «b[i + !(u( 1 ) + J 2 V^Cbu (1) + v( 1 )} + ... ] (4.4) 


(4.1) 

(4.2) 


The second shock condition, (4.2), simply expresses that the tangen- 
tial velocity component is continuous across the shock. Thus, 
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( u (2). u 0);|i + v (2) _ v O) , o (4.5) 

Combining this with (4.4) we find that, to lowest order, 


v< 2 >- v<’> 


Av 

Au 


B + • •• 


(4.6) 


It can readily be shown that the condition (4.2) of continuous tan- 
gential component, to the order considered, is equivalent to the 
requirement that the potential be continuous across the shock. 


Considering first the mid-field region, r = ord (l) , we obtain from 
( 3 . 30 ) and the requirement that the potential is continuous across 
the shock 


A 0 = A 0 q - K A(r Q % v Q ) = O (4.7) 

Now, from (4.6) it follows that, to lowest order, 

A(v q + Bu q ) =0 (4.8) 

across the shock. Therefore, making use of the identity 

1 r 2 2 1 

u v = 7-— (v-Bu)-(v+Bu) 

00 4BL v o o 7 v o o 7 j 

we may rewrite (4.7) as follows 

60 o = - f {aCF 2 ^, 7 !)] - -|(v o ^U Bu o ^^) AJr o } (4.9) 

where 

- ‘5e7 R < v .- Bi o> 


F (? »Tl) 


(4.io) 
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may be termed the "generalized Whitham F-function” as it tends to 
the ordinary Whitham F-function for T| -» oo . The transformations 
(3.28) and ( 3 . 29 ) give 

x + Br = T1 - |g (I) - ?)(v o + Bu q ) - M 2 0 q (4.11) 


x - ** = ? ' h (1 - 0(v o - Bu o ) - M 2 0 q 


(4.12) 


Since both Ax and Ar are zero across the shock these relations 
yield with the aid of (4.8) 



(v ^ 1 ^ + B u 
v o o 



|i?(vi l) + Bu; i )) + M 2 A0 f 


(4.13) 


+ M 2 A0q (4.14) 

From (4.13) it follows that AT] is at most of order u q , i.e. 0(e ), 
so that to lowest order the shock wave transition takes place along 
lines of constant T| . The term involving A0 O is of order A§ times 
u q and is therefore only comparable to the first bracket of (4. 14) for 
the near field r Q = (e)» in which case both are negligible. From (4.13) 
and (4.14) it also follows that Ar /r^^ss = O(e^) so 

that the change in r Q may be neglected compared to r Q itself. It 
therefore follows, that to within the approximation considered (4.9) 
and (4.14) may be simplified to 

M) 

^ F(?,Tl)d| = [(F (2 ^) - (f^) ] (4.13) 



5 ( 2 ) - = KV^B r Q ' (F^ 2 ^ - F^) 


(4.16) 
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The left-hand side of (4.15) follows from the definition (4.10), 
since 





(4.17) 


The relations (4.15) and (4.16) express Whitham ' s "rule of equal areas" 
[8] since by substituting F^ 2 '- F^ 1 ^ from (4.16) into (4.15) one 
obtains 


5 (2) 

\ F(!,T|) d{ = F 0>) (4. 

5 (1) 

Thus, by aid of (4.16) and (4.18) one can use the same graphical 
technique as commonly employed in the first-order far-field solution 
to determine the shock wave location in the mid f ield, provided the 
definition (4. 10) is substituted for the F-function. 


Consider next the near- and the far-field regions. For the near 
field region, r = 0(e), some, additional triple - product terms were 
required both in the expressions for x and r and in the v-velo- 
city component because of the slender-body behaviour near a smooth 
body. However, it is easily demonstrated that for a smooth body 
shocks can only occur for (§/r) = O(s^) in which case the triple- 
product term is negligible. For a body with slope discontinuities, 
shocks will occur also for r = 0(e) ^ = 0( 1 ) , but then the flow 
will behave locally near the corner like a two-dimensional one and 
slender-body triple product terms will be negligible for the deter- 
mination of the shock location. Thus, the mid-field procedure for 
determining the shock location will be valid to the required order 
of accuracy for the near field, as well. For the far field, it was 
demonstrated in the previous section that the mid-field solution 
gave a uniformly valid solution- to lowest order. Thus, .the mid- 
field procedure will give shock locations correct to order unity 

4 

with errors in location of order e and shock strength errors of 

g 

order e . Since such errors are likely to be unimportant in prac- 
tice, the added complication of including the higher-order terms in 
the shock wave location for large distances does not seem worth while. 
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5. SECOND-ORDER THEORY FOR TOE FAR FIELD IN A THREE-DIMENSIONAL 
■ FLOW 

We consider the flow around a simple airplane-like form (wing-body 
combination) located in the half -plane x £ 0 with its (pointed) 
nose at x = r = 0 and the fuselage mainly along the positive 
x-axis. A suitable perturbation parameter e is here the thick- 
ness ratio of the least slender of equivalent bodies of revolution 
in the supersonic area-rule sense. 

With Section 3 in mind, a suitable procedure might seem to be to 
first carry out an expansion in the mid field and then extend it 
to the near- and far-field regions. However, it is found that the 
mid-field analysis employing characteristic variables for the 
general three-dimensional case does not lead to a simple explicit 
solution as in the axisymmetric case. This indicates that the con- 
cept of characteristic variables is not of the same usefulness near 
the body as in the axisymmetric case. Therefore, we will restrict 
ourselves to the study of the far field. Following Section 3, we 
define the far-field domain by 

x - Hr = 0 ( 1 ) 

2 Br = ord e - ^ 

An expansion similar to the one used in the axisymmetric case is 
assumed 

4 i T 8 TT 

u = e U 1 + e U 2 + . . . 

4 8 

v = e V 1 + e V 2 + ... 

4 v 8 „ ' 

w = e W. + e W„ + . . . 

u (5.1) 

x-Br = x=X 0 + e X i + • • • 

4 4 

2 B e r s H = H + c H 1 + ... 

4 ' ' 

9 = 9 q + e 9 1 + . . . 

where U^, V^, W^, X , H^ and 9^ are considered functions of 
order unity of the characteristic variables 5, T| and Q. The ex- 
pressions (5* 1 ) could then be substituted into the equations of 
motion expressed in characteristic form and a sequence of equations 




2k 


generated thereby. However, to obtain the lowest order terms it 
is actually simpler to use (2.l) - ( 2.k ) employing the new variables 
X and H. In these, the irrotationality condition (2.2) becomes 


av ^ h au au 
— . 2E e gH - B ax 


(5.2) 


from which it follows that, to lowest order, 

V, = - B U 1 (5.3) 

Also, (2.3) gives that 


w 1 = 0 


(5.4) 


With the aid of these results, we obtain through substitution of the 
series ( 5 . 1 ) into ( 2 . 1 ) when higher-order terms are omitted 

k„ iZi 

5 1 *X 0 * dH o 


+ 2H U 1 
o 


= 0 


( 5 . 5 ) 


This equation is most easily solved by means of the Monge theory. 
It was shown in [ 10] that the solution can be given in parametric 
form as follows 


U, 


•5^(5, e o ) 


(5.6) 


Xo * 5 - K/H^(5,e o ) 


(5.7) 


By identifying the parameter with § we get a simple and obvious 
relationship between and the Whitham F-f unction. Further, as 

in the axisymmetric case H q must be a function of T] » alone, and 
we may set 


H o = T| v (5.8) 

e o = C (5.9) 
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Next, we proceed by substitution of the series ( 5 . 1 ) into ( 2 . 9 ) - 
(2.14). Using the above results for the lower-order terms we obtain 
the following set of equations to determine the next-order terms: 


3§ 0 


(X e + H,) = 0 


(5.10) 


ax 

an 




aH, 


(5.11) 


_L (m 2 v Tm 2 kT r it \ ( ill > K „ dV 1 

2B 97) (M V 2 + [M -R] bu 2 ) ' 0 2 ( aj / " B H o SH q 3£ 


(5.12) 


ill 2K_L^1 ^o 
an " " b h q ac an 


(5.13) 


(BU 2 . v 2 ) = - ^ 


K 8V . V 1 aH 1 


2B 35 H 35 


(5.14) 


^ 1 BH (M 2 - f) _ V. BH 1 V.. H 1 

BT] ^ BU 2 V 2 ' " H q 3T] B V 1 + H q 31] 


H 


V 2 » H o 1 8H „ 3B 2 2K 3V '- 
+ — -rz— + 


H o an H o an ac ^o 


* 1 \ Uii o s0 1 9V 1 

Vac / an ‘ an ac 


(5.15) 


^2 

a? = " h o bc an 


(5. 16) 


These equations differ from those in the axisymmetric case only by 
the underlined terms, and their integration will proceed in very 
much the same manner as in Appendix D. The final results read 


H i = - x G + n,(n. c) 


(5.17) 


m 2 k - 
M - T n> 


Xl = 


1 ~ B 


K 


\ V (s\ n, £ )ds’ + -^2 ( K + 2 m2 - ) V 1 2h o H o " 

J 2 b M 


.S , 2 


JB V 1 H 1- K ^^2 + *1 


v p- 3 v. 

^ a£ 


( 5 . 18 ) 
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4k 

1 B- bC 


+ ^U.C) 


(5.19) 


BU 2 + V 2 = 


- 1 v i 2 + iH v i(e‘.n.e)<at 


( 5 . 20 ) 


r 


„2 K 
M + 7T 


H, V, 


V 0 = - 


b” v i 2 - -5T 1 * W- 5 v/s'.D.Odt' 

n n o 


bK(1,C) - ,3 V 




2K /_1\ J_ C 
B \9C / ~ H J 


a v 


'.Tl.C)d? ' 


° o a C 


(5.21) 


Wo 


p av. 


H 0 S ac ^ ^»£)d? 


(5.22) 


The underlined terms are the new ones which appear due to departure 
from axisymnetry . A comparison with the axisymmetric solution 
therefore suggests that we may write the far-field solution for the 
three-dimensional case in the following simplified form 


(a) 

u - 

,a u 

2K U ) 

(5.23) 


2 


(a) 

v v J + 

2 “(^) 

(5.24) 

w 

o 


(5.25) 

x<*> + 

K r dw 


o o 

2 a; 

(5.26) 

r<*> - 

aw 

K o 

2B r o d£ 

(5.27) 

t 

S U 


c- ■ 

o 

ac 

(5.26) 
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where superscript 
symmetric case and 
of 


(a) means the same expression as in the axi- 

u , v and w are obtained from the solution 
o o o 


2 S 2 0 o i 

" B 2 + r dr 

. 5x o o 

o 


1 * 2 *o 

„ 2 + 2 2 
or r d C 

o o 


= 0 


(5.29) 


The terms in the expressions for u and v containing derivatives 
with respect to the angular variable £ can be eliminated by making 
the further transformation 


du 

9 o = £ - 2K tf (5.30) 


which leads to the final expressions 


u 

= 

(a) 



(5-31 ) 

V 

= 

(a) 

v' ' 



(5.32) 

w 

X 

= 

w 

o 

(a) 

Kr o 9w o pT , K 2 / Su o> 2 

+ 2 ae + 2BK r 0 \d 9 ) 

O O 


(5.33) 

(5.34) 

r 

= 

(a) 

Kr dw ,du 2 

-2B°30 - 2K r o( 3 6 ) 

o x o 


(5.35) 

0 

= 

e - 

o 

du 

2K ^ 


(5.36) 

Here , 

u o’ 

v and w are now solutions 

o o 

of ( 5 . 29 ) with 

£ replaced 

by 0 

o 

• 

In considering the shocks in the 

far field one 

notices 


that to lowest order the solution has the same functional behaviour 
as the axisymmetric one, except that 0 Q appears as additional para- 
meter in the F-function. One can therefore use the same method as 
in the axisymmetric case to determine the shock wave location to 
first order. 
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6. EVALUATION OF THE F-FUNCTION FROM WIND-TUNNEL EXPERIMENTS 

A new method to determine the F-function from wind tunnel measure- 
ments, was proposed in [ 14] . In it, flow deflections are measured 
instead of pressures ns in previously used methods. Thereby , con- 
siderably improved accuracy is made possible, particularly at 
high Mach numbers, since flow deflections are much easier to mea- 
sure accurately than pressures. The relationship between flow 
deflections and the F-function becomes simple only in the far field, 
so that the flow deflections must be measured sufficiently far away 
for the flow to have approached its far field behaviour, but yet 
close enough where the deflections still reasonably large to allow 
good measurement accuracy. Inspection of the linear solution for 
axisymmetric flow reveals that its far-field behaviour is attained 
when Br is large. For slender configurations the far field 
character may be dominant already for values of Br greater than 
3, or sometimes even less. This makes the new method particularly 
useful for high Mach numbers for which the far field may be approach- 
ed already for r = ord(l). 


The new experimental technique utilizes measurements along a cyl- 
indrical surface r = const, of flow deflections a and t in 
the radial and azimuthal directions, respectively. From figure 6.1 
one obtains 


tg a 


v 

1 + u 


tg T 


w 

1 + u 


1 


> 


( 6 . 1 ) 


e 

r 



Figure 6 . 1 
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In order to make direct use of the results from the previous section 

to relate these to the F-function we shall tentatively assume that 
/ -4 s 

Br = ord (e J so that these results are directly applicable. Through 

expansion for small a and t and making use of the far-field 

approximation for u, v and w in terms of u , ’V and w q one 

can derive expressions for u , v , w in terms of n and t 

o o o 

(see Appendix E). These, in turn, can be used in the expressions 
relating the velocities to the F-function producing the following 
relations : 


p (?' 9 o> - /?<*.•> -isfEl 

I — ' 0 

+ ( B + J if cr2 ( x > r ) ( 6 *2) 


where 


§ = x - Br + 2Kr a (x,9 ) 

,2 


/. ,2 K\ x 
(M +3-) 


B 


^ CT(x',e)dx' - 


(|f^) + |r CT 2 (x,e){B 2 -1-|- 


- [K + ^ M 2 - In 2Br } 
M J 


(6.3) 


and 


0 _ Q _ 2K d<j(x;9 ) 

o B 90 


(6.4) 


Here we have omitted the additional terms that arise for a body of 
revolution with slope discontinuities. 


The error in F would be, formally, of order e 10 if, as assumed 
in the derivation, r = ord(e - ^). However, when used for moderately 
large values of Br, the error will of course be larger and prim- 
arily governed by the neglected terms in the linear solutions. The 
terms which are quadratic in o (underlined in the formulas) are 
therefore small compared to the others as they essentially represent 
higher-order cumulative effects for large r , and could therefore 
usually be ignored. The expressions presented in [ 14] are then 


recovered . 



30 . 

If; the F-f unction evaluated in this way is plotted as a function 
of § for constant 0, then the presence of shocks implies that 
it will be undefined for certain ^-intervals . However, because the 
F-function is. used only in the study of the far field, one- can - as 
a consequence of Whitham ' s • "rule of equal areas" (see Fig. 6.2) -• 
overcome this difficulty simply by defining the F-curve in each 
of the "missing" intervals as the straight line connecting the 
end points of the curve. . 



a = b + e 
a + c = b + d 


j- ■ C' + e = d . 


Figure 6.2 
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7. CONCLUSIONS 

By the use of characteristic variables a second-order solution valid 
in the mid field r = ord (l ) of an axisymmetric flow was derived 
which was found to be identical to the one obtained before [ 13 ] by 
a different method. The mid-field solution was then extended to the 
near field, r = 0(e), basically by adding terms equivalent to the 
triple product term in Van Dyke's [7] solution. Also, correction 
terms were added to the expressions for the dependent variables so 
as to make the transformation of the body contour continuous when 
the body has slope discontinuities. The mid-field solution is valid 
to first order in the far field, r = ord (e^);by adding further terms 
which are of higher order, in the near field it was made valid also 
to second order in the far field. Thus, a second-order solution 
for axisymmetric flow uniformly valid everywhere has been produced. 
The solution is expressed directly in terms of the linear solution 
and may be considered an extension to second order of Whitham's [8] 
first-order rule. An extension of Whitham's method for construct- 
ing the shock wave was also presented and shown to give the shock 
strength to second order in the mid field. As the new solution 
also allows the presence of body surface slope discontinuities, it 
should be of practical value for the calculation of the supersonic 
non-lifting flow around slender bodies of revolution even in such 
problems in which the sonic boom is not of primary interest. 

For the general three-dimensional case, introduction of character- 
istic coordinates does not seem to lead to any substantial simplifi- 
cations for the near and mid fields. For the far field, a second 
order solution was obtained through the new procedure. The charac- 
teristic coordinates used employed the downstream and upstream 
Mach conoids as surfaces of constant § and T| , respectively, and 
perturbed meridian planes Q = const, intersecting the downstream 
Mach conoids along bicharacteristics. Ihat the bicharacteristics 
perhaps are not so fundamental in the three-dimensional case is 
indicated by the fact that the simplest asymptotic far-field repre- 
sentation, i.e., the one resembling most closely that for axisym- 
metric flow, is obtained by instead using for the third angular 
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coordinate 0 , where 0=9 + A0 and where A9 turns out to 

o o u 

be half the angular displacement of the bicharacteristic along the 
downstream Mach conoid through that point [.13]. 

The asymptotic three-dimensional far-field solution can be put to 
practical use in at least two obvious ways. First, its simpie 
structure indicates how good approximate solutions can be generated, 
and some suggestions were given earlier [ 13 ] how one could use it 
for such a purpose. Secondly, it makes possible the generation of 
simple formulas relating flow deflections in the far field to be 
Whitham F-function which can be used in a new experimental proce- 
dure to determine the sonic boom strength from wind tunnel measure- 
ments. Application of this method to a simple body of revolution 
reported in [ l4] shows that it is capable of yielding accurate 


results 
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CHARACTERISTICS IN THE AXISYMMETRIC . CASE 

For axisymmetric flow the system (2 . T.) — (2 .k ) becomes 


(«* _-<■*) g * ^ - = 2 > 

3 V 3U 

3x - 3r 1 ' ' 


3U \ 
3 r ‘ 



J 


(A.1) 


When studying -such a system-, it would certainly be an advantage if 
x-derivatives and r-derivatives did nof'appear in the same equation. 
Theref ore, we should be interested in the possibility of finding for 
example a new set of independent variables T| such that the 

system could be brought into a form in which one of the equations 
contained derivatives' only with respect to £ and the. other, one . 

a ■* * * 'fr • ’ * ' ' , . / ^ 

only with respect to T); Now it is shown in some text-books (see, 
e.g. , [4] , p. 101 or [6], p. 433) that iri every domain where 
the flow is supersonic (i.e. q > c) the system (A.l) can be brought 
into such a form, namely 


3U v 3 V 

3? Cot( * + a) + 

3U ./ x 3 V 
— cot(#- a) + — 


c 2 V 3r 
V 2 - c 2 r ^ 
c 2 V 3r 
^-c 2 r 9T 1 


\ 

► 


where the transformation 


U,Tl)- (x,r) 


is given by 


(A. 2) 


If - ft »<*<*- «> 

i 

3x 3 r . / „ \ 

■ TH = ai) cot( * + a) 

and c, $ and a are to be found from 


(A. 3) 


tfu 2 a V 2 ) 



V 

1 


u tg $ 
sin a 



\ 

> 


j 


(A.4) 



3 6 

The system (A. 1 ) is then said to have been put into charac t er is t ic 
form with ^ and T) as characteristic coordinates . The curves 
^ = const, and T| = const, in the xr- plane are known as character - 
istics (downstream resp. upstream characteristics). (See [4], 
pp. 62-7 2, for example). 

If we introduce the disturbance velocities u and v and make 
use of (A. 4), then (A. 3) becomes ( 2 . 3 )- (2.6) and (A. 2) becomes 
(2.7) - (2.8). 
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APPENDIX B 

TRANSFORMATION OF THE THREE-DIMENSIONAL FLOW PROBLEM INTO CHARAC- 
TERISTIC FORM 

With Appendix A in mind it seems reasonable to try to introduce 
new independent variables £ j T| and Q such that the system 
(2.1)- (2.4) is brought into a form in which at least some of the 
equations would contain derivatives with respect to only two of 
these variables. The new variables § , T| and Q will be called 
characteristic coordinates. 
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We notice that 


and 


a . . 
Ji 


9x. 


a . . 

ij 


q ± 


i Further , u and 


'\( u 1 


u 


2’ 


^ 3 ) 


3_ 

3 5 


are general notations for the 


and ( 


3 v • 


3 

3 u. 


3 u. 


-) , respectively. 


'/the matrices 


row vectors 
We also introduce 




It is easily seen from the above and Bernoulli's equation that a i j » 
A^, h and H are functions of Q and r only. 


By means of these notations we can now write, the system (2.1)- (2.4) 

*. ’ . • • f • 

and the equation (B.2) in a more compact form, namely 


I 

i= 1 


3Q 

A. = H 

x 3x^ 


and 


I 


a . . 


3 2 j> 


i, j=1 


ij 3x,3x. 
1 J 


= h 


(B.3) 


(B.4) 


respectively 



With a coordinate transformation 
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y ± = y i ( x 1 > x 2 < x q) 


3' . 


i = 1,2,3 


(B.3) and (B.4) — expressed in the new variables y^, y ^ and 
appear as 


3 


l 


i= 1 


3Q 

ay ± 


= H 


and 


I 


b. . 


a 2 <& 


i, j=1 


ij Sy^yj 


= k 


respectively , 


where 


B 


- 1 r^ A - 

1 L 3x . j 

j=i J 


b . . 
Ji 


^ ^y ± Qy-j 

k- - ~ Zj ^1 dx , dx . 


ij 


k,l=1 


and 




i , j >k, 1=1 


\k 1 


B y k 9x l 
ij ax ± 3xj 9y k 1 


a . . 


Suppose now that we could find a function ij/(x) satisfying the 
order equation 


I 


i. j=1 


a aj[_ djr_ 

ij 9x i Sxj 


= 0 


Then by choosing, for example, y^ = f(x) we would have b^ = 


y 3 

(B.5) 


(B.6) 


(B.7) 


(B.8) 


first 


(B.9) 


0 



2 

from (B.8) which means absence of from (B.6). Consequently, 

dy-| 

if Q were known on only one surface = const., it would not 

be possible to determine ty~ uniquely from (B.5) on that surface. 


Such a surface is known as a characteristic surface (see [6] , 
Chapt. VX, § l). Thus the surfaces ty(x) = const. — where | is 
a solution of (B. 9 ) — form a family of characteristic surfaces. 


In our special case we easily find that for arbitrary real numbers 
P-|» P 2 3X1(3 P 3 

3 

V a ij Pi Pj s Q + (x,p) • n_(x,p) (B.io) 

i»j=1 


(B.11) 

(Here U, V, W, c and hence also a„ are considered as given . 
functions of x ) . 

Thus (B.9) can be written 


where 


ft ± (x,p) 3 Up 1 + Vp 2 + ¥ — ± c 



n ± (x, ^-) =0 (B. 12) 

9x 

I 

From (B.11) we conclude that (B.12) has a solution only if M T > 1 . 
Now we consider only the case with > 1 everywhere in the flow 

field. Evidently (B. 9 ) is a partial differential equation of the 
first older. Within the general theory of such equations (Monge's 
theory) (see [4] , Ch. II or [ 5 ], Ch. 2 or [ 6 ], Ch. II ) the 
concept of characteristics plays a very important role. The charac- 
teristics belonging to (B. 9 ) are called the bicharacteristics of 
(B.4) (see [6], Ch. VI , § l). The equation for an arbitrary bi- 

.characteristic follows from Monge's theory and reads 

dx = t!!± ( - a± ) 

dt gp ’ 9x 

where t is a curve parameter and ijf satisfies (B.12). 


(B.13) 
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It. is easy to verity that \|r is constant along such a curve, thus 
every bicharacteristic lies entirely on a characteristic surface 
ijr = const. One may observe that, in case \f. 0, equation 

(B.13) is equivalent to 


dx 2 9n± { ~ M. \ 

d — = ( x ’ at } 


dx. 


2 
d.flfct 


1 = ( 

' 4 -*- dp, ' 


d t 


C it 


dx 


l ) 


ij; = const. 




(b.i4) 


Consider for a moment the case of. undisturbed flow,: i.e. 
which is a trivial solution of (2.1)- (2.4). 


Q = U 


By putting 


*1 

^2 


= x - B r 

= x + B r 

= 0 


one easily finds .that 


- ) = o 


n_( 


sn_ _ 9yi 

and .. ( x , — — ) = 0 




dx ■ 


i = 1 ,2 . 


Thus (see (B. 12) and (B.l4) the surfaces y^ = const, and y 2 = const 
are two families of characteristic surfaces, the first of which is 
intersected along bicharacteristics by the surfaces y^ = const. 

y ^ = const, is known as a downstream Mach cone , 
y^ = cons t , as an upstream Mach cone and 


const, is a meridian plane 



k2 

Now we return to the general supersonic case, only with the restric- 
tion ..that the disturbance, caused by the body, is considered to be 
small. In the choice of new independent variables we can then take 
over some ideas from the special case above. 


Zl is chosen such that : 




/ 3 1 \ 

y 1 = const, is a characteristic surface (0_(x, and a 

disturbed downstream Mach cone (x - Br w const.). 


^2 

r 2 


is chos en such that : 


By. 


2 ' = 0) and 

disturbed upstream Mach cone (x + Br ss const.). 


y 2 = const, is a characteristic surface (ft _ (x, -■ = — ) 


*3 


is chosen such that : 


y ^ = const, intersects y 1 = const, along a bicharacteristic (see 

(B.I 5 ) below) and is a disturbed meridian plane (9 ss const.). 

From (B.14) we get the condition for y^ = const, to intersect 
y 1 = const, along a bicharacteristic, namely 


dy^ = dy 1 = 0 => 


or 


dx 2 



dx^ 



3Q_ 

(x, 

) 

3x 

9 ft 

(x. 

a 7l 

ap 2 

3P^ 

9x 

ax 2 



Bx^ 



Sy 2 



sy 2 



3ft_ 

Bp 2 

(x, 

) 

dx 

3 ft_ 
3 p 3 

(x,. 

By 1 

dx 


= 0 


= 0 


(B.15) 


Thus, for the determination of the transformation (x^x^x^)-* (y^y^y^) 


we have the following equations : 
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“ (B. 16) 
(B.17) 


(B. 18) 


( (B. 18) is obtained from (B. 15 ) by means of (B. 17 )). 


We now return to the system (B. 5 ), which is equivalent to the dyna- 
mical equations (2 . 1 ) - (2.4) expressed iri the new variables.' Because 
y 1 = const, and y 2 = const, are characteristic surfaces we get from' 
(B.8) and (B. 9 ) that b^ = b 22 = 0. But this means,, as mentioned 
before, that dQ/3y^ (i = 1,2) cannot be uniquely determined from 
(B.5) on a surface y^ = const, on which Q is known. 

We therefore try to rewrite (B. 5 ) such that one of these equations 
does not contain dQ/dy 1 and another one does not contain 9Q/by 2 * 

By putting 



one obtains from (B. 7 ) and (B.9) that 
X^ B ± = 0 , i = 1,2. 

Thus multiplication of (B. 5 ) from the left by and X^^ gives 

the scalar equations 



X 4 (2) *i) & - 


i=1,3 


X 


(B.19) 


(B.20) 


i=2 ,3 

respectively. 

To these equations we add the last two from (B.3) expressed in the 
new variables 


v ( \ 

L \ dx„ dy, / q 1 

k=1 J 

V ( _JL_ V 

L \ dx dy, ) q 2 

k=1 3 k 


- X ( 


3x 1 9y k 3 


X 

k=1 * 


(B.2 1 ) 


(B.22) 


(B. 19 ) - (B.22) are then equivalent to the original system (2.1)- (2.4) 
expressed in the new variables y 1 , y 2 and y^. 

We note that it is not necessary to set up a potential equation (B.2) 
in order to define the concepts of characteristic surface and bicharac- 
teristic as has been done above. We could just as well work directly 
with the system (2 . 1 ) — (2.4) and obtain the same results, but we 
follow here the usual approach taken in standard textbooks dealing 
with second order scalar equations instead of with systems. 

Now we put 

t - r, ' 

n = y 2 ■ 

C “ y 3 , 

Noting that 


9y ± 3(y r y 2 >y 3 ) 


dXj ^(x^XgjX^, 


Sy i-1 


9y i-1 


9y i + 1 


9y i + 1 


i , j mod . 3 
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(B. 16) — (B.22) can be written 
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(B.24) 
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In these equations 

_L 

c f, V-1„2 ro 2 2 2-, \ 2 

a ( 1 - - T g' M[2u + u + v + W ] 1 

00 

which is easily obtained from Bernoulli's equation. ‘ 
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APPENDIX C 


EXPANSION OP THE LINEARIZED AXISYMMETRIC SOLUTION FOR LARGE r 

The solution of ( 3 . 25 ) can be written (see for example [8], 
p. 307) ' 


x - B r. 
o o 


u (x ,r ) 
o v o’ o' 


f(t) dt 


o '/(x -t ) 2 - B 2 r 2 




x -Br 
o o 


1 (x o -t) f(t) dt 

r (x ,r ) = — \ — 

°°° r ° ^ Jf 7*2 2 

o . v (x -t) - B r 


(C.1) 


(C.2) 


.. 


with f given by the integral equation 


x - BR (x ) 
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( x o _t ) dt 

i(x n -t) 2 - B ¥(x o ) 


= R (x ) 10 (x ) , x s 0 (C.3) 

O' O'' o' o' ’ o ' ' 


where u) q (x o ) is of order e (see 3.47)* 


Now (C.3) gives for f, approximately, 


f(t) « ^ [R Q (t) «J 0 (t )3 


Hence 


f = O(e^) 
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In the following we will use the variables 5 and r Q instead of 

x and r 
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For 0 £ | < 2 B r Q we can expand ^ 1 + ^ ^ 

and integrate term by term. 


in a series 


Thus , observing that 


t 1 t 
n- 1 n 


$ f (t) (?_t) * dt = $ dt i S dt 2 • • • $ [S dt l dt ” 

o o n 


o o 


we obtain for 0 ^ £ < 2Br 


S' (- 1 ) 11 r ( 2 n-l)il 1 F n^) 

U ° ' ' 4 n! L 2 n J ( 2 B r ) n+ * 

n — n * n ' 


> (C.5) 


. v = - B 
o 


V 2n + 1 (-1 ) n r ( 2 n-l); ’ 1 F n^^ 

L 2n.1 n! L 2 n J ( r y 

n — r* ' n ' 


n+2 - 


where 


3 

F*(?) = ^ F n-l( ? ’) d?< ( n = 1 » 2 , 3 ,...) 


►.(«) ■ ^*( 5 ) - 


From (C.k) and (C.6) one easily deduces that 

F n^ = [( 2 n+ 1 )i J f ( H g ? )] ‘ ? 2 

^ ✓ 

0(e 2 ) 


(C.6) 


(C.7) 


for some in the interval [0,1], 
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APPENDIX D 


THE FAR FIELD IN THE AXISYMMEIRIC CASE 


Assume the expansions 


= e 

U 1 + 

e U 2 + . 

4 

8 

= e 

V 1 + 

• 

+ 

CD 

B r 

= X 

= X + e 

2B r 

s H 

= H + e** 



O 

r k’ V k 

’ *k 

and 


order unity and = ord (1 ) « 


These are substituted into the system ( 2 . 5 ) - (2.8) which leads to 
the following set of equations: 


To the lowest order we get from ( 2 . 5 ) - (2.8) 

( 


3H 

c 


= 0 


dy 1 9H „ 
o 1 of /\,2 


+ TT 


an 2 an 


[(M 2 - K) u, ♦ V,] ■ > 0 


3 

95 


V 1 3H 
1 ■ o 


(B U 1 + V r ) + — ^ = 0 


H q 9? 


Mj < BU 1 - v l> - if 


V. 3H 
1 o 


o^ 


= 0 


(D.1) 


CD. 2) 


(D.3) 


(d.4) 


Equation (D. l) implies that 


H c = h o (n) (D.5) 



Substitution of (D.l) into (D.3) gives 


B U 1 + V 1 =0 


Further, by introducing U 1 from (D.6) into (D.4) we arrive at the 


equat ion 


. av 1 aH o 

■ 2 W " % aF = 0 


which has the solution 


V 1 = B 


^(%) 




(D.7) 


By means of (D.6) and (d. 7) we find from (D.2). 


Fl- -U 


x 0 = * 0 (5) -k/h^^(?) 


in which, in order for the solution to 'take - the correct, ..foiw for 7 - . 

zero disturbances, we choose £ 0 (5) = 5 * 

Using (D.5) - (D.8) we obtain from (2.5) - (2.8) to ; the next order":;; 


If. ^o + H,)- = 0 


(D.9) 


3x 1 _ _K_ m 2 i 40 - 2 _K_ ~"1: . .... 

%T " h Ji L 3 “ ‘ m 2J a-n ,v i . " 2 b. an-i- vr , 




' 25 arf (m2 v 2 + ^ ^ B > 


a 3V/ V. 3H 1 

af (B u 2 + v 2 ) = - 2B T|- - ir.aT 

3 1 3H M 2 - | „ 

V 1 3H, V, Hl 3H o • V 2 3H o 

+ H o 3T) ‘ H 2 + H 0 

o 
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From (D.9) we have 


H 1 = " X o + h l( 1 l) 


Using (D.7) - (D.9) we get from (D.ll) 


B U 2 + V 2 = 


I v i + 1T S V 1<5M) d ^' 

O *•' 


Substitution of U 2 from (D. 14) into (D. 12) leads by aid of 
to the equation 


o o o 

3 H ^ 

5^5if J V i (5 ' ,!,)d5 ' 


i a «. ’f. *1 .2 

H ST) B V 1 
o 


from which we find, using ’(D.7), 


V 2 * 


^f„2 H 1 V 1 3 1 


B 


v i 2 - 2jr lt ifrS v t (?, ’ 11)d5 ' 

o o ° 




By means of (D.7), (D. 14) and (D. 15) we get from (D. 10) 


*1 = 


m 2 k 5 

M - if 


B 


- C V (g',T])d§' + “2 (K + | M 2 - ^|) V 2 H In 
«J 2B M 


- 2B V 1 H 1 - K ^o^2 + 


Putting 

2 ^ \ 6 


e ^( 5 ) + c°^ 2 (0 ■- F(5) 


e 4 ^ = X(S) + f F 2 (§)[m 2 - 2 + |+4(K + jM 2 -p)ln 


(D.13) 

(D. i4) 

(D.7) 


(D.15) 

*0 " 

(D. 16) 

] 



53 


and remembering that 


H q = H - e 4 H 1 + 0(e 8 ) 

we finally obtain from (D* 6 ) — (D. 8 ) and (D. 14) — (D.16) 


u = 


7^, + ^-372 \ F (5 ' )d? ' + (M 2 - f ) £ f^ L - + o(e 12 ) (D. I?) 

S 2Br 4*(2Br) J/2 ^ 2 2B ^ 


v = - B u + 


B 


(2Br) 


372 ^ F(g')d^« - BK L f - ^- + 0(e 12 ) (D. 18) 


x - Br = 5 - K \f 2Br F(^) + 


M 2 -! 5 

/2Br 


^ F(?')d5' 


+ § (K + | M 2 - ^|) [P(§)] In 2Br + | (M 2 - 2 + f) [F(?)] 2 + 
M 
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APPENDIX E ■' ' 

DERIVATION OF u q , v q , w q , 0 ‘ AND t q /t AS FUNCTIONS OF 

a, t, AND r IN THE FAR-FIELD DOMAIN (r = ord (e~ 4 )) 

By means of (5*3l)— (5-33) s (5«35)> (6.1) and the formulas 

Bu o + v 0 + = °( e12) ’ 

• , ' . - ° • . _ • • ; 

8 

0 = 0 o - K r o u o v o + 0(e ) and 

-Bu + o(e 8 ) = : -Bu q + dr(e®') = 

= v Q + 0(e 8 ) = v = 0(e 4 ) , 

which can be deduced from ( 5 . 17 )- (5.22), we obtain > 

v = ( 1+u)tg ct = ^1 - cr + P(e ) ’ = ( 1 _ Jg ■[( 1 - o + 0(e 12 ')j-^ o + o(e 12 ) . 

V v = (1 - -g a) a + 0(e 12 ) . 


v + 0(e ) /; M + K \ w' 1 2 \ 

' = p 1 = ( 1 + 5 — v v + 0(.e ) 

° 1 +(M 2 + K) u q V B J 

( „ B 2 + K N 12 

• v o = (J + — b cj) a + 0(c ) . -• 

X X 

i = ^ u dx ' = - | ^ v dx 1 + o(e 8 ) . 


V 0 = - -g ^ a dx' + 0(e 8 ) 


0 o = 0 + K r o u q v q + 0(e 8 ) = 0 - v 2 + 0(e 8 ) 


o o 

X 


V 0 O = - ^ ^ a dx' - | r a 2 + 0( e 8 ) . 



B 2Br 


* ' IT ‘ IS ♦ °<‘' 2 > ■ 


b 2 *| 


u = (1 + M 2 u o ) U Q :+ 0(e 12 ) = u o + ^ V 2 + 0(e 12 ) 


i K 
- 1+ 2 


= (1 + u) tg T = T + 0(e 2 ) 


w = t + 0(e ) • 


w = w + 0(e ) 

o ' ' 


w o = T + °(e ) . 


g- = 1 - | V •+ 0(e 8 ) 

1 - K u q + 0(c 8 ) B 


.. r o . K 8 , 

■ — = 1 - g o + 0(e ) • 


u = " I ( 


1+ — x 
o 


+ o(c 12 ) 


= ( 1 - b a ) ° + °( c 2 ) 


n , 12 •. 
= t + 0(e ) 


1 ^ ct dx' + 0(e 8 ) 


0 
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and 


u = 
o 


-a 


i + 


b 2 + | 

B 


) CT + — l — dX ' + °(® ^ ) 

' ZB T ^ 


\ 


( , B 2 + K \ . 12 v 

v o = V + ~B“ CT ) CT + 0(6 } 


w q = r + 0(e 12 ) 


V (E.2) 


0 = 
^ o 


I 5 CT 6x1 - 


K 2 8, 

■g r ct + 0(6 ) 


r 

o 

r 


- g a + 0(e ) 


y 
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